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1 Introduction 



The nature of the dimensionality of the manifolds appearing in euclidean quantum 
gravity ensembles has attracted considerable attention recently IQ @, d, H H 0- 
There are, in principle, many ways of defining the dimension; on smooth regular 
manifolds we expect every definition to yield the same value (for example 2 if the 
system is "two-dimensional"). However the euclidean quantum gravity ensembles 
contain many members which are very far from smooth, often showing some kind of 
fractal structure, and different definitions of dimension probe different global aspects 
of the geometry. Thus they can (and sometimes do) yield different numerical values. 

Although these questions are of very general applicability by far the most com- 
plete studies of the notion of dimension have been conducted for two-dimensional 
euclidean quantum gravity (here "two-dimensional" means that every manifold in 
the ensemble is of dimension two in terms of the standard mathematical definition of 
a manifold) . The ensemble is defined at fixed volume (which we denote the canonical 
ensemble, CE for short) by the partition function 



Z(V;X) 



J Dg5^J fty/g-Vjexpi-S^ftoX)) (1) 



where the functional integral runs over all physically inequivalent metrics g with the 
volume constrained by the delta-function to be V. The functional S e ff(g; A) is the 
effective action obtained for the metric after integrating out all matter fields and A 
represents any parameters (for example the central charge c) describing the matter 
field theory. Expectation values are given by 

UcE = z^j Dg(-)6^j d 2 ^g-V^e^{-S eff { g] \)) (2) 

with the understanding that it only makes sense to calculate the expectation values 
of reparametrization invariant quantities. We can also define a grand canonical 
ensemble (or GCE) partition function 

/■oo 

Z(fi; A) = / e-^ v Z(V] A) dV. (3) 
Jo 

In practice nearly all the available results about the dimensionality of the manifolds 



in these ensembles have been obtained in the discretized formulation |8|, |9], [I(J . The 
functional integral over metrics of a given volume is replaced by the sum over all 
triangulations (or, more generally, graphs), Q, with the number of vertices, Ng, fixed 
so that the CE partition function is 

Z(N; X) = J2 W e exp (-&//(£; A)) (4) 

Q 
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and expectation values given by 

^ CE = Z{N-X) Z^^(") exp (~^/( g;A ))- 

To calculate the GCE partition function the integral over V in ([|) is replaced by a 
sum over N. 

Two different notions of dimension have been considered in detail. The first is 
the "Hausdorff dimension", dh, which is defined by considering the volume of space, 
dr, contained in a shell of geodesic radius R and thickness dR; provided that R is 
much bigger than the short distance cut-off and much smaller than the characteristic 
linear geodesic size (which is V 1 ^ for some exponent v) we expect that 

(dr) CE ~ R dh ~ l dR. (6) 

For c = an analytic calculation has been done and it is found that dh = 4 ||. For 
other central charges c > the situation is still unclear but numerical results indicate 
that d h = 4 for unitary c < 1 theories || and there is some analytic support for this 
0| . For c > 1 it is believed that the ensemble of universes collapses to a branched- 
polymer (BP) like phase; for a pure BP ensemble analytic calculation shows that 
dh = 2 |l]J which is in very good agreement with numerical simulation for c = 5 
(which seems to be the value of c above which all these results change very little). 

The second notion of dimension which has been investigated is the "spectral 
dimension", ds- This is defined for a given smooth metric g on a d-dimensional 
manifold by the behaviour of the coincidence limit of the heat kernel K g (t\ £, £') for 
which it is known that 

1 oo 

K g (t^,0 = ^J2 fa r(0 ( 7 ) 
r=0 



where the functions a r (£) are reparametrization invariant, see e.g. [12|. Thus in the 
quantum gravity ensemble we can define ds by 

for t small enough that the random walks generated by the Laplacian do not probe 
the finite size structure of the manifolds. In the discretized formulation we consider 
a random walk on a graph Q with iV sites 0; the random walker moves from one 
lattice point to one of its neighbours at each time step with equal probabilities 
corresponding to diffusion governed by the Laplacian of the graph. Letting the 
probability that the walker has returned to the starting point % after t steps be 
Pg(i; t) we can define ds by 



p N {t) = (N-'j^Pg^tyj 



~ t d S /2- ( 9 ) 

CE 
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This behaviour is expected to hold provided <^t N 1 ^ for some exponent A so 
that discretization and finite size effects are avoided. Numerical simulations show 
that ds ~ 2 for c = and that as c increases d$ decreases 0. It was conjectured 
in 0] that du = 2g?s for all values of c. This would imply that if the large c phase 
is a BP phase then it has ds = 1. The purpose of this paper is to compute ds 
analytically for a pure branched polymer phase in order to check this conjecture 
and provide some definite result with which to compare the numerical simulations. 

This paper is organized as follows. In section 2 we explain how ds which is 
defined in the CE can be calculated in the GCE. In section 3 we define the branched 
polymer ensemble and find relations satisfied by the return probabilities. In section 
4 we do the GCE sum and show that ds = §■ Section 5 contains the proofs of various 
results used in section 4; some technicalities are relegated to the appendix. Finally 
in section 6 we discuss our result and its connection to other spectral dimension 
problems (in particular on percolation clusters) that have been considered in the 
past. 



2 Spectral dimension in the CE and the GCE 

Our calculation of the spectral dimension for the BP ensemble (which we will define 
carefully in section 3) depends crucially upon the relationship between the CE and 
the GCE. This is most easily understood in the following way (we assume from now 
on that we are working with the discretized formulation). Consider first the CE; 
at very large times t ^> N A , where A is a suitable exponent, the walk will have 
uniform probability of being at any of the N lattice sites (essentially this is just the 
contribution of the zero mode of the Laplacian) so 

lim P N (t) = 1 (10) 
At much smaller times we expect the behaviour (|9|) so that 

A suitable interpolating function between these two limits is 

^ = jf + tsa^{-w) (12) 

where a is a constant. From this we can compute the generating function 

oo 

P N (y) = £ y'PMt) *>x— y + (1 _ y + x--y-^ (13) 
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where we have assumed that ds < 2 and only written the dominating terms in 
Pn{v) as y f 1 and N — > oo. Note the appearance of the simple pole at y = 1 
as a consequence of the zero mode and that the rest of Pn(v) is analytic at y — 1 
for finite iV. Defining P' N (y) to be Pn{v) with the contribution of the zero mode 
subtracted (which we will find is straightforward in our actual calculation) we see 
that 

so that, as we might expect, the non-analyticity associated with dg appears only in 
the thermodynamic limit. Note that we also have that at large N 

P' N (1) ~ N A{1 - ds/2 \ (15) 

This result depends only on there being a cut-off in the t~ ds / 2 behaviour at t ~ iV A 
and not on its particular form. 

In the GCE the natural quantity to compute is 

oo 

g t=o 

where the sum runs over all graphs Q in the ensemble, Ng is the number of points 
in Q and, as before, Pg(i;t) is the return probability of a walker on Q returning to 
site i after t steps (in general one would also sum over starting points i but in the 
case of branched polymers it is sufficient to consider walks starting at the root and 
summation over i is not necessary). The function V is related to the CE quantity 
Pn(u) through 

V(z,y) = J2^Z(N;X)P N (y). (17) 

N 

We expect that for large N the CE partition function behaves as 

Z{N; A) ~ N Jstr ~ 2 ZQ N (18) 

where zq is a constant and ^ str < \- It follows that the coefficient of the pole term is 
finite at criticality (ie at z j zq where graphs of arbitrary size contribute). Dropping 
the pole term to get 

V'(z,y) = J2^Z(N;X)P' N (y) (19) 

JV 

and approximating the sum over iV by an integral we find, using (fOp, 



( \ 

i-y 



(20) 
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where the prefactor exponent (3 is given by 



(3=1- lst r + A ( % - 1 ) (21) 



and the function <3> by 



Mr) / ae .^ . - dx. (22) 



Note that $(t>) is analytic for v > with $(0) = 0(1). The result (f2~0|) describes 
correctly the singularity of V(z,y) at z = z provided that the prefactor diverges 
as z | 2 ; if it does not then we consider instead a higher derivative of V'(z, y) with 
respect to z. We draw three lessons from (pD|) : 

1. It is not the case that 

In fact, for z < Zq, V'(z,y) is an analytic function on a neighbourhood of 
y = l. 

2. The exponent A is a gap exponent; that is the n + 1st derivative of V(z, y) 
with respect to y is more divergent than the nth derivative as z f by a 
factor (1 — z/,2 )~ A . 

3. If we can compute the prefactor exponent, (3, and the gap exponent A in the 
GCE we can find ds using (|2~T|). 



This discussion is based on the interpolating ansatz (|T2| ) and the reader might worry 
that it is not completely general. In fact, although the ansatz is useful for pedagog- 
ical purposes, we shall show in section 4.2 that it is not necessary. 

3 Branched Polymers 
3.1 The Polymer Ensemble 

In this work we will restrict our attention to polymers which have vertices of order 
one or of order three and in which one vertex of order one, called the root, is 
distinguished. This is the simplest case of the generic rooted branched polymer 
ensemble ||11|| . The simplest polymer, Bo, consists of one link joining the root, 



labelled "0", to one other vertex of order one which we label "1" as shown in fig.l. 
The ensemble of all polymers, £>, is generated by the statements: 
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Figure 1: The elementary polymer Bo- 

1. The only element of B consisting of one link is Bo- 

2. Any two polymers, Bi and B 2 may be combined to generate a new, and larger, 
polymer B 3 = Bi U B 2 as shown in fig. 2. 




Figure 2: The combination of polymers B x and B 2 to form B 3 . 



It is clear that any polymer B e B, B ^ B , has a unique decomposition into two 
constituent polymers B' and B" by cutting just to the right of vertex "1" as shown 
in fig. 3. 




Figure 3: The unique decomposition of polymer B to its constituents B' and B". 
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One characteristic of an element B of B is the number of order one vertices (not 
counting the root) which we will call the size and denote Nr; clearly we have 



iV BlUB2 =N Bl +N B2 . (24) 
Letting be the number of elements of B with a given size, N, we find 

AT-l 

QjV — ^m^n-m (25) 

M=l 

with Qi = 1. This is the recurrence relation for the Catalan numbers, Tjv, so 

(2iV-2)! 

^ = T ^ = iV!(iV-l)! - (26) 
For large N, the numbers T/v have the asymptotic behaviour 

Tjv ~ -3-. (27) 

JM 2 

Note that the number of vertices of order three is N — 1, and hence the total number 
of vertices excluding the root is 2N — 1, so N is a sensible measure of the volume 
of the system. We therefore expect on general grounds that Vt N ~ N lstr ~ 2 exp fj, N; 
comparing with (^) we recover the well known result for branched polymers that 
•jstr = |. The generating function of the Catalan numbers will appear repeatedly in 
the rest of this paper and is given by 



oo 

T{z) = * N Tn = r(l - Vl^Tz). (28) 



2 

N=l 

Note that this definition of B naturally connects ensemble members of different sizes 
suggesting that the Grand Canonical Ensemble (GCE), which contains polymers of 
all sizes, may be easier to calculate with than the Canonical Ensemble (CE) which 
contains polymers of fixed size; indeed (^) is nothing but the GCE partition function 
for B with the action iVe- 



3.2 Random Walks on Polymers 

Since B is defined with a vertex of special status, the root, it is convenient to consider 
walks starting and ending at the root rather than summing over all vertices. This 
does not affect the value of ds because in fact half of all vertices have coordination 
number one like the root and are therefore equivalent as starting points. First 
consider an arbitrary polymer B ^ B . At t = a walker sets out from the root; at 
t — 1, after taking one step, he/she is necessarily at the first vertex; for the next step 
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one of the three links attached to the first vertex is chosen with uniform probability 
and the walk proceeds down that link. The process continues in this way with the 
step taken chosen with uniform probability from the possible steps. Now let P B {t) 
be the probability that the walker returns to the root for the first time after t steps, 
and let -Pb(^) be the probability that, after t steps, the walker is at the root. Note 
that both Pn(t) and -Pb(^) are zero for odd t and that P^(0) = because at least two 
steps are needed to return to the root for the first time. Clearly walks contributing 
to -Pb(^) may be back at the origin for the first, second, third .... occasion so we 
have 

n 

P B (t)=5 t , + P£(t)+ £ P ^ tl )Pi(t 2 ) + ...+ £ Y[pi( ti ) + .... (29) 

h+t 2 =t ti+t 2 +...+t n =t i=i 

Introduce the generating functions 

oo 

p B (y) = J2y t/2p ^ t ) ( 3 °) 



t=Q 



and 



t=0 



We substitute into (pOQ and then use ([H]) to obtain 

P B (y) = l + Pi(y) + (Pi(y)) 2 + ■ ■ ■ 

= i=w < 32 > 

Let us now consider a polymer A ^ Bo which is made by joining two polymers 
B and C, cf. fig. 2. We want to relate the first return probability on A to the first 
return probabilities on its constituents B and C. A typical walk leaves the root and 
then goes out and back many times from the first vertex along B and C before 
finally returning to the root; an example is shown in fig.4. Taking all these walks 
into account we find that 

Pi® = \<>t,2 + Q) 2 (Ph{t ~ 2) + P&t - 2)) 

+ (T) 3 E ( P B(tl)+P^)(pl(t 2 )+Ph(t 2 )) + ... 

^ ' t!+t 2 =t-2 

•••+ 3 E \{{Pi{U) + P l c {U)) + . . . . (33) 

h+...+t n =t-2 i=l 
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Figure 4: Example of a walk contributing to the first return probability on polymer 
A = B U C. 

Inserting ( |33"D into (pHD we get 

PKv) = ^ + (|) 2 ( P B^) + P c(2/)) + (|) 3 (^b(2/) + ^c(2/)) 2 + - • • 

•••+(£f +1 (Pb(v) + P^(y)T + ■ ■ ■ 

!l (34) 



3-(P£(y) + P£(y))" 
It is convenient to introduce a new quantity h^(y) denned by 



hx{y) = (l - • (35) 



We find from fl32|) that the return probability generating function on A is then given 

by 

Pa(v) = ^— ( 36 ) 

i - y h A {y) 

and the recurrence relation for the first return probabilities (|34|) becomes 

hA[y) = n-(i-y)(My) + M!/))- (37) 

9 



Note that for the elementary polymer, Bo, we have 



h Bo = 1. (38) 

We can, at least in principle, generate Pb{v) f° r every B 6 B from the recurrence 
(|37|) and the initial condition (j3S|). 

4 Random Walks in the GCE 

4.1 The generating function 

We define the GCE generating function 

Q(z,y) = J2z NB PB(y) (39) 
which, because of fl36|), becomes 

e(«,y) = T ^-E 2HVb (Vy))" 1 - ( 4 o) 

First we will show that Q(z, y) has a simple pole at y = 1 with finite residue 
provided z < \. From (]37f) we see easily by induction that the functions /ib are non- 
decreasing. It therefore suffices to show that the sum in ([40"D with /ib(?/) replaced 
by /ib(1) converges. From (|37D we have 



/lBU0(l) = l + Vl)+Ml)- ( 41 ) 



From this, the initial condition (PB|), and the relation between the sizes (23), it 
follows that 

/i A (l) = 2iV A - 1. (42) 

Note that this quantity takes the same value for all polymers of a given size. Using 
(f|^) we have 

(i-y)Q(z,y)\ y =t-J2w^l (43) 

N=l 

which converges for z < \ on account of (P7|). 

As we discussed in section 2, the simple pole in Q(z, y) at y — 1 arises from walks 
that are much longer than the size of the polymer; we expect to find these walks at 
an arbitrary point in the polymer with essentially uniform probability which in this 
case, according to (|42|), is (2A^ A — l) -1 - The information about ds is contained in 



10 



the remaining (ie non-pole) part of Q(z,y) so we define a new function 

d 



Q(z,y) 



dy 



((l-y)Q(z,y)) 



Now define 



Quiz) 



BeB 



d 



-2 dh B (y) 
dy 



(44) 



BeB 



d 



y=l 
n+l 



dyj h B (y) 



(45) 



y=l 



Assuming for the moment that all these objects actually exist, then formally 

G(*.v) = £^^Q»(*)- (46) 

n=0 

We show in the appendix that this series is asymptotic but that it is not Borel 
summable. Nonetheless, we will argue in section 5 that the leading terms can be re- 
summed to an integral representation which is analytic in the region of interest. The 
absence of Borel summability means that such a function is not unique. However 
it can only differ from the true result by terms which have essential singularities at 
y = 1 (and vanish there) and we do not expect these to affect the spectral dimension; 
so long as there is any power law scaling as y — > 1 any vanishing essential singularities 
are irrelevant. 
Now let 



h 



(n) 



d 

dy 



hv{y) 



y=l 



Then 



d 
dy 



n+l 



n+l 



= > ^77^ > (n + l;ai, 

y=l r =l \' l B I T(n+l,r) 



n+l 



0-n+l 



(0> 

13 / 



(47) 



(4* 



i=l 



where (to; a\ . . . a m )' is the multinomial coefficient [13 

to! 



(to; ai . . . a m )' 



(l!) a i . . . (m!) am ai! ...a m \ 



and the summation region ^(n+l, r) consists of all non-negative integers cti, 
such that 

n+l n+l 

yka k = n + l, y^ j a k = r. 

k=l k=l 



(49) 

O-ri+l 

(50) 
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Inserting fl48[) into fl45j) and reordering the sums gives 

n+1 / rr n+1 r/) (i VA 

r=l ^(n+l.r) VBgB J / 



fl-{ni,n 2> ...,n p ) ^ = J^Z* 13 fj/i^ , p, 71* > 0. (52) 



Let us define 

Bee j=i 

We will prove in section 5 that 

^(ni,n 2 ,..,n p )^j ^ ^(m ,n 2 ,...,n p ) M _ 4^|-p-| ELi n « (53) 

where the symbol ~ means the most singular part as z f | and where y4( n i> n 2,...,n P ) 
is a constant. Using we note that 



^(ai®l,a 2 ®2,...,a n+1 (g>n+l)^ __ (<lz — \\ | IL=1 (^4^)°" J 

V cfe ' \Be/3 (^B ) r+1 / 

where by g ® m we mean that the number m appears q times in the list. Using 
fl5"3| ) and integrating r + 1 times we find the leading singular behaviour of Q n {z) for 
n ^ to be 

a. w . (1 _ 4zr i-g (z^ nfe) 

^— f 2 T+i I(^n + r + l) 

r=l v z ' 

(n + l-a 1 ,...a n+1 )A {a ^ 1 ' a2(S,2 '---' a "+ im+1 \ (55) 



x 

^"(n+l,r) 



When n = we obtain 



&(z) = -~log(l-4z). (56) 



Substituting (|55|) and fl56|) in ([46]) we see that the most singular part of Q(z, y) as 
z | | is the sum of a logarithmic piece and a function of (1 — y)(l — 4z)~i. For the 
rest of the analysis it is more convenient to work with which takes the form 

9S 1 M (57) 



dz l-Az \(l-Az) 



for some function From now on in this paper we will be working with this most 
singular (as z \ j) part of Q unless otherwise stated. Note that from (|57| ) we deduce 
that the gap exponent A = §. 
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4.2 Extracting ds from the GCE 

The behaviour of Q(z, 1) = Qq{z) given by (|56|) can be used to place an upper 
bound on ds which we will need later. From (f2~7|), (|39"D , ( f44[) we have 

Q(^,l) = ^iV-i(4^ J P;(l). (58) 

Together with ([56]) this implies that -Pjv(l) ~ ATs at large AT. But from (|Uj) this 
can only happen if ds < 2. 

As we discussed in section 2, the GCE generating function is related to the CE 
generating function by 



dz u \ u 



poo 

I dx x l3tr ~ l $(1 — y, x) exp(— xu) 
Jo 



(59) 



where, for convenience, we have introduced u = 1 — Az and a continuous variable x 
equivalent at integer values to N, the size of the system. Since we are looking for 
the asymptotic behaviour at large A" we can approximate the discrete sum over A" 
by a Laplace Transform in x; the difference will be sub-leading. The extra factor 
of x appears in compared to (IE) because we are looking at ^ rather than 
Q; this makes life easier because it removes the possibility of a trivial divergence at 
small x. From (|59"D it follows that 

$(1 -y, x ) = x l ~ lstr (f)({l -y)x A ) (60) 

for a suitable function <\>. To show this introduce new variables u' — u(l — y)^ 1 ^ 
and x' = x(l — y)^ A into to get 



^ ( Ja) = C 1 - l/)~ (7 ^- 1)/A dx' x'^-^(l - y, (1 _^ )1/A ) exp(-xV) (61) 
but by setting y = and it = it' in (pS|) we also have 

— $ ( ] = / dxx 7s "-" 1 $(l,x)exp(-xn'). (62) 
u' V M / Jo 

Similar identities are valid for the derivatives of $ w.r.t. z so 

$(1, x) = (1 - y )-(7^-D/A$(i _ ^ _^_ ) (63) 

and the result flSOf) follows. Now it is easy to see that $(1 — x) must be a bounded 
quantity in x > for fixed y < 1; we have 

1 Tjv 00 / 1 \ 

•(l-».W) = 5r-EEWfl«-MfTT) < 64 » 

JV j=l t=o v 7 
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where the subtraction is to remove the pole term. Now Pi(t) < 1 and -Pi(O) = 1 so 
we get 

It follows from fl65]) that lim I _ >00 $(l — y, x) must be finite for y < 1 and this is 
compatible with (|60D only if <p(s) behaves like, 

(j)(s) ™ s -(i-7-r)/A ( 66 ) 

and so we conclude that 

lim $(1 - x) ~ (1 - ? / )-( 1 -t^)/ a . (67) 

Substituting in the values of ■y str = \ and A = | for the branched polymer ensemble 
B we obtain 

lim $(l-y,x) ~ (68) 

— 2 

This is the generating function liniAr^oo P' N {y) so we find that lmiAr-^oo P/v(t) ~ t~s 
and hence that ds = § (see ([IS]) and (|H"D). 

4.3 Sub-leading terms 

There are two sources of subleading behaviour. The first is the corrections to the 
asymptotic form (B6|) which must be like 



s 



) ™ g -d-7*)/A ( ! + £^f*^ e>0 (69) 



in order for (^) to be satisfied. The correction term is more singular as y | 1, 
and so has smaller ds, but is suppressed by N e in the thermodynamic limit. This 
is expected because we know that there are members of B (for example the linear 
polymers which are basically infinite chains with finite sized outgrowths) which have 
ds = 1 but are thermodynamically insignificant. 

The second source of subleading behaviour is the corrections to (^) which, as 
we will show in section 5, take the form -u 5+ 5~P~f SLi n » with 5 = |, 1, | . . .. The 
consequence of this is that (1571) is modifed to 



9Q 1 _5 / l-V \ , 1 £ / 1-3/ 



^En^M^ (70) 



1-42 \jl-4z)t,/ ^(l-4z)^ \(l-4^) 

where the functions $j(l — AT) have the same upper bound as $(1 — y,N), see 
(|65|). Repeating the manipulations of section 4.2, and assuming that the z'th term 
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in (ffOD makes a contribution which survives as N — > oo we find by the arguments 
of the previous sub-section that it has a spectral dimension 

4 = 2 + ^(25-1) > 2. (71) 

This shows that it is not possible that the leading term, i.e. the $ term, makes no 
contribution as N — > oo and that the lower bound ( |55] ) is fulfilled by one of the 
subleading terms; this would lead to ds > 2 for B in violation of the bound ds < 2 
obtained from the behaviour at y = 1. 

5 Proof of the Principal Result 



In this section we will prove the result fl53|) and derive constraints on the coefficients 
A^ nx '— np '. It is convenient for what follows to adopt a more concise notation. Let S 
denote a list of non- negative integers ni, . . . n Ps and let 



We adopt the convention that the empty list corresponds to ps = 0. In this notation 
we will prove that leading singular behaviour as z f j is given by 

PS 

BeB i=l 

~ A^u^-ps-hs^ u =l-4z, (73) 

where A^ is a constant. We will refer to the power of u in ( |T3"D as the actual degree 
(of singularity) of H^(z). If ps = we have the special case 

H0(z)=T{z)*± (74) 

We note that, having established ( |73"D for a given set of integers S, it is a corollary 
of (|73D that 

PS 

H^ s \z) = $>^(/4 0) ) fK l) , m>0, 

BeB i=l 

* A^ V ^ + VS + q -^ u^s- q -hs. (75) 
This follows from fl4~2] ) since 



BeB i=i 
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HH v,(i 



1 d_ 
2 



y ^(ni,n.2,...,n p ) u |-ps-|ns 



(76) 



and the result follows. 

We start by computing the derivatives for polymer A = B U C in terms of those 
for its constituents B and C; differentiating n times and setting y — 1 we obtain 



h 



(n) 
A 



4 n) + >4 n) 



n! / ci 



^ n — r! yc^/ 



(/ iB (2/)) r (l + / iB (2/)) + (B^c) 



+E(l)(^d/)) fc (^(i/)r fc + 



A:=i 



j/=i 



Consider first the case p = 1, ni = 1. We have 

fcg> = + fcg> + 4 0) (1 + ^ 0) ) + hP(l + &g>) + 2hg>hg\ 
Summing over all A = B U C (and remembering that = 0) we find that 

AeB 



HW(z){l - 2T{z)) = 1(1 - Az)-I + 0((1 - Az)- 1 ] 



(77) 



(75 



E^K ) +4 0) (i+^ 0) ))E^ c + ( B - c ) 

= 2T(z) (# (1) (^) + H m (z) + H^(z)) + 2H {0 \z)H {0 \z). (79) 
This can be rearranged to give, using (|73|) and (|75|) , 



(80) 



Note that the most singular term on the r.h.s. of (|80|) comes from the H^°'°\z) term 
in (|79D . Finally, substituting for T(z) from (p8|), we find 



51) 



in accordance with (1731). 
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We can now proceed to the general case by induction. The strategy is always 
the same as the one we have just used. We take the definition 

PS 

H^(z) = J2^ A U hi A ) (82) 

AeB i=l 

and replace by its expression fl77p in terms of its constituent polymers B and 
C. By doing this for sets of integers S chosen in a particular order we can ensure 
that the resulting expression on the r.h.s. of fl82f ) contains only H^(z) which are 
already known. To see this we can write (|77| ) as 

h<£> = /4 n) + /4 n) + T(n - 1) (83) 

where by T(n — 1) we denote all contributions which contain derivatives of no higher 
order than n — 1. Let S = (rii,n 2 , . . . ,np s ) where each n* < n — 1. Then consider 

H (l®n,S)^ = J2z NB J2 zNC ( h B ) + h C )+T ( n -Vy 

BeB CeB 

PS 

x J] (C> + + T{n t - 1)) . (84) 
-i=i 

After multiplying out and doing the sums over B and C the generic term on the 
r.h.s. of (jS4]) is of the form 

H {Sl \z)H {s *\z) (85) 

where S\ and 5*2 are lists of integers with the property that no member is higher than 
n. We will denote by S*_ a generic list of numbers derived from S by removing one 
element and by S + a generic list of numbers derived from S by removing one element 
rik and replacing it by {mi, ...} where the mi add up to n^. We can categorize the 
possible pairs of lists S\ and S 2 that appear in (|^) by keeping track of the number 
of factors of T which appear on multiplying out. 

1. No T factors. Apart from the two terms which reproduce the l.h.s. multiplied 
by T(z) these give S\ and 5*2 each with either a) fewer than q elements with 
the value n and all other elements less than n (since they are drawn from S) 
or b) q elements with the value n but fewer than p$ remaining elements less 
than n drawn from S. 

2. One or more T{n — 1) factors and any number (up to p$) of T(rii — 1) factors. 
In this case each of Si and S2 contain n fewer than q times with all other 
members being less than n. 
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3. One or more T(rii — 1) factors and no T(n — 1) factors. In this case Si and S2 
contain n up to q times, at most ps — 1 members drawn from S and the rest 
are drawn from the integers < n — 2. 

From la and 2 we see that H^ q ~ l ® n ' s \z) must be calculated before H( q ® n > s \z) 
and from lb that H^ qtSn,s -\z) must be calculated before H^ n,s \z). Finally, from 
3 we see that H^ q ® n ' s+ \z) must be computed before H^ q ® n ' s \z). These criteria 
are automatically satisfied by lists of numbers ordered in the following way (below 
S a > S b means that is calculated before H^^). Let the largest number 

appearing in the list S a be n a and denote the number of appearances by q a . The 
ordering is defined as follows: 

1. If n a > n b then S a > S b . 

2. If n a = n b , and q a > q b then S a > S b . 

3. If q a = q b then replace n a and n b by the next largest numbers appearing in 
S a and S b respectively and apply l.,2., and 3. recursively until the issue is 
decided (if there are no more numbers in a list then the corresponding n is set 
to zero). 

Of course to calculate as far as a given S it is not necessary to know the result for all 
preceding lists in the ordering defined above because ([84]) only requires knowledge 
of preceding lists with at most ris + 1 elements. 

Now we need to keep track of the degree of singularity of the terms ( jS5D appearing 
on the r.h.s. of (|84j) . Define the naive degree of as 

doih^) = -1 - \n, (86) 
and the naive degree of a product of such terms 

doiflht^^itdoih^). (87) 

i=l i=\ 

Then, if (|73| ) is already established for the lists Si and 5*2, the actual degree is given 
by 

d = 1 + d , (Case 1) (88) 

unless it happens that the product consists entirely of and no h^^ (or vice 
versa) in which case 

d = I + d , (Case 2) (89) 
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(this is because, anomalously, ~ ~ and not tta). Now compare two terms of naive 
degree d' Q and d' ' respectively which satisfy rf' < d' f — |; then the corresponding actual 
degrees tf and <i" satisfy ff < c?" because 

max(d') = dJ, + l<t^ + - = min(cf ). (90) 

This result can be used to discard terms in ( jS2|) that can only make sub-leading 
contributions to H^(z); only terms whose naive degree is at most \ greater than 
the lowest naive degree appearing in an expression can possibly contribute to the 
leading divergence at z f 4. We will call this result the "naive degree criterion". 

To keep account of the naive degree in subsequent formulae, we will adopt the 
notation that quantities in square brackets [ ] immediately after a term denote its 
value of do. Anotating (|77|) in this way 

hf = ^[-1-In] + (B^C) 



+ E^TH (dy) {(hMY^r-ln} + {hMY^r-l^n] + (b-c) 

+ Il( r k )( h c(y)) k (h B (y)y 



^- fc ri r+ i _3 n l 



+ E( r t 1 )(^(y)) fc (^(y)) 

k=l ^ ' 



■ (91) 

y=l 



Using the naive degree criterion we see immediately that there is only one term in 
the braces which can ever contribute to the leading divergence which can therefore 
be computed from the expression 



PS / 

H( s \z) ~ J2 zNB J2 zNc Ii[ h 



PS 

K) 



L> l- 1 -!"! 

L'y-_B C6B i=l X 

Ui — 1 

[-|-§ n ] + ( B ->c) ). (92) 

2/=l 



When we multiply out the r.h.s. of ( |9"2"D the naive degree criterion again tells us that 

.1 _ 3 
2 2 



only contributions containing at most one factor of the term with do = — \ — |n 



need be retained leaving 

H {s \z)(l-2T(z)) ~ H^\z)H^\z 

SiUS 2 =5 

PS "fc-1 



+2T(z) YsY. nk \ ) H^- l - m ' m ' s ^ (z) (93) 

k=l m=0 \ m / 
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where the primed sum indicates that the empty list is not included, and S\rik is the 
list S with the number n& removed. Both terms on the r.h.s. of ( |93D have actual 
degree 1 — ps — \n$ and so the result ( [75| ) follows for H^ s \z). 

Having established the leading behaviour of it is straightforward to cate- 
gorize the possible sub-leading behaviour. From fl93|) we note that one contribution 
to the first sub-leading term of arises when on the r.h.s. we take the first sub- 
leading term of H^ Sl ' and the leading term of (or vice versa). Assuming that 
the first sub-leading term of takes the form u a ~ bps ~ cns and equating powers we 
get 

which immediately tells us that b = 1, c = |. From (^l|) and the result (|80|) for H^ l > 
we also see that the first sub-leading term is just vh less singular than the leading 
one. Repeating this argument for lower divergences we conclude that a general 
sub-leading contribution to must take the form 

M 5+i- Ps -|n S; <y= .... (95) 

We can also use (f^|) to obtain an expression relating the coefficients 

PS n k -l 



SiUS 2 =S k=l m=0 ^ ' 



■m,m,S\n k ) 



(96) 



Although these recursion relations cannot be solved exactly we can use them to 
obtain estimates of the coefficients A^ s >. First we show that all A^ are positive. 
Assume that the A^ are computed in the sequence described above. From (ff5|) 



we see that if A^ is positive for some S then so is A^® '^ and it follows that if 
A& > then also A^+) > 0. Since A® > (from (gT|)) it follows that all 
are positive. It is convenient to define a new coefficient B^ s \ which is also always 
positive, by 

ps 

A(s) = B ^l[n t \ (97) 
i=i 

and substitute this in (|9~6|); all the factorials cancel leaving 

Ps n-k—l 

B (S) _ B(Si) B (S2) + B (n k -l-m,m,S\n k )_ (ggj 

SiUS 2 =5 k=l m=0 

We can obtain a lower bound on by introducing B^ which satisfies = 
BO, B^ = B^°\ and 

PS n k -l 

B (S) _ B (n k -l-m,m,S\n k ) ^-gg-j 

k=l m=0 
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which is (|98[ ) with the quadratic term dropped. Now (^) is solved by = 
b(n s ,p s ) where 



b(n,p) = nb(n — l,p + l). (100) 
Iterating this equation we find that 

b(n,p) = n\b(0,p + n). (101) 
But b{0, p + n) = B((H-»)«o) which we can compute exactly from (ff^), so we obtain 

6(n,p) = 2^- a n! r(n ^ ) "^ ) . (102) 

Now note from ( [75] ) and (ITJ2) that 

£ (5 ' g ® 0) > &(n S ,p g + g) (1Q3) 
b(n s ,ps) 



This ensures that those terms on the r.h.s. of (|9q) and fl99|) which arise when m = 
or rifc = 1 are larger for the true coefficients B^ than they are for the B^ s \ It now 
follows that, since ( PD| ) is just (P§D without the quadratic term, 

B^>b(n s ,p s ). (104) 

Finally note that, as a consistency check, for p = 1, ni = 1, this lower bound 
evaluates to give = ~ which is the exact result (pTj) . 

To find an upper bound on B^ we first need the subsidiary result that 

Q = Yl T ( an Si +PSi ~ |) r («^5 2 +Ps 2 ~ \) 
Si us 2 =s 

< 4T(l)T(an s + p s -±). (105) 

This can be proved by replacing the Gamma functions with their integral represen- 
tation and using the inequality w + s@ < (t + s)" for t, s > and (3 > 1. Now 
define 

S^f ) = 5 r(a n 5 + p 5 - I) r Qns+ps (106) 
where B , a, and r are constants. By fll05|) we have 

J2' B^B^ < 4T(l)B B { + ] (107) 

Sius 2 =s 
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and 

PS Wfc— -1 

^ B M-,™,%) = ( ^ n . j 5o r(an s + p 5 - i + (1 - a)y- a r an ^ 



k=l m=0 




(an s + ps - \ + (1 - a) - (ps - \ + 1 - a)) 



a 

xT{an s +ps-\ + 0-- a))r ans+ps 

< B ° rl a T(an s + Ps - i + (2 - «))r awg+ps . (108) 
a 

If we choose a < 2 then this is generically bigger than .ETj. whereas if we choose 
a > 2 it is generically smaller. Combining ( |105| ) and ( 108 ) we find that for a = 2 

PS n k -l , v 

^ B^B™ + B^- 1 -™'™'^ < ( 4T(i)S + — ) flf } . (109) 

SiU5 2 =5 fc=l m=0 ^ ' 

Now note from (f75|) and ( [106|) that 

-wsr<^r' < 110 > 



This ensures that those terms on the r.h.s. of (|98| ) and which arise when m = or 
rifc = 1 are smaller for the true coefficients B^ than they are for the B+\ The 
numbers B^ provide consistent upper bounds on B^ if we can choose Bo and r 
such that the r.h.s. of ( |109| ) is less than B+^ and such that B^ < B+ . It is easy 
to see that this can be accomplished by choosing r = 2 and 

B = ( - + e] -4tt (HI) 
where e is small and positive so we conclude that 

B (s) < B T(2n s +p s - \)2 2n ^ . (112) 

Both the upper ( |112j ) and the lower ( |102| ) bounds on B^ behave as T(2ns+ P s~ 
i) with an exponential prefactor. With coefficients of this form the Taylor series for 
Q(z,y) (|46|) is re-summable which we will demonstrate for the upper bound B+\ 
From PB| ) and ( |5oT ) we have 

Q(z, y) = ~ log(l - Az) + G(z, y). (113) 

Replacing the coefficients B^ by their upper bounds B+^ and using ( p7[ ) and ( |106| ) 
G(z, y) becomes 

G + (zy) = -T( V ~ l Y"-V(-ir +1 Hr(|n) 2 2 ^> 
G+( ^ 2^^(i_4,)i; n»2J 1} nin + r + l) 2 

x V (n + 1)!r(2(T2 + 1)+r -^ (H4) 
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Using the integral representations of the Beta and Gamma functions and interchang- 
ing the order of the summations gives 

„ . d f 1 dt f°° e~ s ds f .„ . 5 ^ / 3 9n „ 

G ^ v) = W. tJL — {(i-o-»E(-- 2 r 

+I ^ ( - ir ^(-^r £ (»*) 



where 



r=2 m>r JF(n+l,r) 



u> = ^ — J-,. (116) 
(l-4z)§ 



The remaining sums can be done leaving the integral representation 

„ a r l dt r e- s d S 

G+{Z ' V) = dw~ L 7 



o Jo 



J (i _ - s - 2t " (e^-l-C)} (117) 

L 1 — wt^s 2 ) 



where 



C = -*(l-t) Wth l • (118) 

1 — wt^S 2 

We see that, as claimed in section 4.1, G + is a function of w. In the region of 
interest, 0<y<l, 0<z<j, wis always negative and it is straightforward to 
check that the integrals are finite, not only for G + but also for all derivatives of 
G + with respect to w. Thus G + is analytic in this region; it is a function whose 
asymptotic series about y = 1 reproduces the series for G with the coefficients 

(S) 

replaced by B\ . It is not at all trivial that this is the case; indeed if the coefficients 
were to have the behaviour T{2ris + /3ps — 5) with [3 > 1 then the resummation 
would lead to a divergent integral. 



6 Discussion 

The value ds — 3 can be compared with the numerical simulations of ]7[ and in 
particular fig. 10 of that paper. The data for central charge c = 5 is in good agreement 
with the analytic result for B. The c = 5 data also yields a value of j s t r close to \ and 
a very good fit to du = 2 which is the analytic result for the Hausdorff dimension of 
B. Taken altogether this is very strong evidence that at c = 5 the quantum gravity 
ensemble is truly branched-polymer like. On the other hand our result clearly rules 
out the conjecture that dh = 2d$ for a unitary quantum gravity ensemble (this 
conjecture is also inconsistent with very high precision numerical results for c = —2 
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which is a non- unitary theory ||). For values of c closer to one the numerical results 
for ds are some distance from ~ just as those for dh and ^ stT are some distance from 
the generic branched polymer values. This is a familiar phenomenon in these models 
which seem, at least from the numerical point of view on finite size systems, to have 
complicated cross-over behaviour from c < 1 properties to generic large c properties. 
This can arise from competing contributions which obscure the asymptotic large 
volume behaviour and has been discussed at length in [13]. Recently simulations 
have been performed to determine the spectral dimension of what is believed to be a 
branched polymer phase in four-dimensional euclidean quantum gravity; the results 
are in excellent agreement with the expectation that ds = 3 [pl| • 



The calculation presented in this paper shows that ds = § for the simplest 
example of a generic branched polymer ensemble. By tracing the origin of the gap 
exponent A it is straightforward to see that the same structure will emerge for any 
other generic ensemble. The crucial point is the factor (1 — 2T{z)) = (1 — Az)^ 
appearing on the l.h.s. of ([53); the power | just arises from the structure of the 
GCE partition function and is simply 1 — •y str where, generically, 7 str = |. For 
the non-generic branched polymer ensembles which can be produced by allowing 
vertices of arbitrary order and weighting them in particular ways [IH] it is possible 
to generate different values of 7 str and we expect that these ensembles will in general 
have ds 7^ §. 

Finally we would like to discuss the connection between our calculation and the 
spectral dimension on percolation clusters [JIT]. These authors originally conjectured 
that the spectral dimension on percolation clusters at criticality is | independent of 
the embedding dimension D. High precision numerical calculations have long since 



shown that this is not correct for small D (see [[18| for a review). However for infinite 
(or sufficiently large) D this result is believed to be exact and there are a number 
of scaling arguments and approximate calculations for it |TjJ. Now consider the set 
of (bond) percolation clusters (PC) on a Cayley tree with vertices of order three 
and with the constraint that the bond from the root to the first vertex is occupied 
(fig. 5). This set is very similar to the B ensemble except that the percolation clusters 
contain two-point as well as three-point vertices. It is straightforward to map the 
clusters onto a branched polymer ensemble which contains graphs with two-point 
and three-point vertices with appropriate weights. The only subtlety arises in the 
weighting of the two-point vertex. Every time a two-point vertex is present in 
a branched polymer there are two percolation clusters corresponding to it (fig. 6). 
Thus the branched polymer ensemble which generates the percolation clusters with 
the correct weights has GCE partition function Z satisfying (see fig. 7) 
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Figure 5: Percolation cluster on a Cayley tree. 



Z = z{l + a 2 Z + a 3 Z 2 ) (119) 

where a 3 = 1 but a 2 = 2 and z is the bond occupation probability in PC language. 
This modified branched polymer ensemble still has •y s t r = \ and hence, according 
to our discussion above, ds = §; in fact it is the simplest extension of £>. Because 
of our interest in quantum gravity we defined d$ in terms of probabilities averaged 
over all clusters of size N rather than in terms of a particular "typical" cluster. 
The calculation picks out the leading non-analyticity as y — > 1 after taking the 
thermodynamic limit; we have 



lim $(1 - y, N) ~ — — — = - — L— (120) 

This shows that in the thermodynamic limit a finite proportion of the Tjy clusters 
of size have ds = |. As we discussed in section 4.3 clusters with dg smaller than 
| are present but the number of such clusters is suppressed by a power of relative 
to Tjv and so the probability that they appear as an infinite percolation cluster is 
zero. 
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PC BP 

Figure 6: Mapping of percolation cluster two-point vertex to branched polymer. 




Figure 7: The branched polymer equation that generates percolation clusters. 
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Appendix 

In this appendix we prove that the GCE generating function (|39|) with the pole term 
at y = 1 removed has an asymptotic expansion about y = 1. 

The return probability -Pb(^) on an individual polymer B can be expressed in 
terms of a transfer matrix Tb = 1 — -D B where is the Laplacian on B. Z} B has 
eigenvalues {A B , i = 1 . . . iV B }; one eigenvalue, A B B , is zero while the rest satisfy 
1 > A' B > 0. We have 

PB(t) = J^a B (l-\ B Y (121) 
i=l 

where the constants a B depend on the eigenvectors of D-q and satisfy < a B < 1. 
Substituting ( |121|) into (|39"D but omitting the zero eigenvalue (which generates the 



pole at y — 1) we obtain 



Q\z,y) = J2 zNB T,y t E 4(1 -A B r- (122) 



BeB i=0 
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Q! is absolutely convergent for 0<y<l, 0<z<j. Using the binomial expansion 
we get 

SQ'(z,y) = Q'(z,y)-Q'(z,l) 

oo t / j_\ Nb — 1 

= E * Nb EE L ) {y ~ 1)n E «b(i - Ay*. (i 23) 

BeB i=0 n=l ^ ' i=l 

The term (y — l) n can be replaced using the integral representation 



x n 



1 f°° 

-- I v^e-dv (124) 



o 



so that, by re- arranging and re-summing over n and t, we get 

POO 

5Q>(z,y)= (3(v)e~^dv (125) 



o 



where 

m = "E^ B E^B^^exp (-v 1 —^) . (126) 
Bee i=i V B / 

We can now use ( |125| ) to develop a formal series in powers of (1 — y) by repeatedly 

integrating by parts which gives 

Q'(z, y) = Q'(z, 1) + (1 - y)P(0) + (1 - y) 2 (3'(0) + . . . + (1 - y) n ^ (0) + i? n (127) 

where 



/■oo 

i? n = (l-l/r / f3 (n \v)e~^dv. 
Jo 



'1281 



To compute the nth derivative of f3{v) from (|126 ) is trivial and so is the remaining 



v integral in ( |128[ ) (remember that y < 1). We find the remainder term 

^ = (i - y)"(-i) n+1 E ^ B E a B ; b) i¥ — ( 129 ) 

BeB i=l — + 

and the coefficient 

iV B -l 

/3(n)( ) = (-1)^ £ ^ £ _ ^n+l^-n-^ (13()) 

Be/3 i=l 

We see immediately that 

R n = (l-y) n+ \-l) n+l r n (131) 

where 

r n < (3 {n) (0). (132) 
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It follows that, provided the j3^ n \0) exist, the series will be at worst asymptotic. 

Letting A™ n be the smallest non-zero eigenvalue and using the bound a l B < 1 we 
find that 

\(3 {n) (0)\ < J2 zNBN B^B in )~ n ~ 2 - (133) 
Bee 

Now A™ n ~ N B 5 where the largest value taken by 5 is 5 = 2 (the value for the linear 
polymer). This leaves us with an upper bound on f3 n (0) of the form 

r(2n + |)(l-4^)- 2 "- 9 / 2 . (134) 

Of course this is a substantial over-estimate compared to the detailed calculation 
presented in the body of this paper but it tells us two things. Firstly the series (|127|) 
is asymptotic. Secondly, since the linear polymer is only exponentially suppressed 
in the ensemble, /3 n (0) must grow much faster than T(n) and hence the series cannot 
be Borel summable. 
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